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In two dimensions the non-interacting density of states displays a Van Hove singularity (VHS) 
which introduces an intrinsic electron-hole asymmetry, absent in three dimensions. We show that 
due to this VHS the strong-coupling analysis of tunneling spectra in high-T^ superconductors must 
be reconsidered. Based on a microscopic model which reproduces the experimental data with great 
accuracy, we elucidate the peculiar role played by the VHS in shaping the tunneling spectra, and 
show that more conventional analyses of strong-coupling effects can lead to severe errors. 
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Scanning tunneling spectroscopy of Bi-based cuprate 
high- Tc superconductors (HTS) shows a d-wave gap and 
a strong dip-hump feature which is nearly always stronger 
for occupied than for empty states .J- It has been pro- 
posed that the dip-hump structure results from the in- 
teraction of electrons with a collective mode^ but the 
dip asymmetry has not received an explanation so far. 
Indeed such a coupling leads to electron-hole symmetric 
spectra in classical superconductorsi^i^i^ The dip-hump 
was also observed by photoemission, but in those ex- 
periments it is not possible to probe the electron-hole 
asymmetry. Sometimes photoemission spectra are even 
symmetrized,^ thus ignoring the relevance of the asym- 
metry seen in tunneling. The fact that in two dimen- 
sions the density of states (DOS) has a prominent van- 
Hove singularity (VHS), unlike in 3D, introduces nat- 
urally an asymmetry and thereby modifies the strong- 
coupling analysis and the corresponding determination 
of the collective mode frequency in an essential way. 

Photoemission experiments have provided a detailed 
account of the band structure in cuprates ii'^>^ In agree- 
ment with early calculations,— the band crossing the 
Fermi level presents a saddle point leading to a logarith- 
mic VHS in the DOS. The scanning tunneling microscope 
(STM) is the ideal tool to look for such singularities, since 
under suitable conditions it probes directly the DOS with 
meV resolutioniiiii^ii^ Up to now, however, there has 
been no report of a direct STM observation of the VHS 
in HTS materials, neither in the normal nor in the super- 
conducting state. Previous interpretations of the missing 
VHS invoked the tunneling matrix element i ^^i^^ in planar 
junctions, it is indeed believed that the DOS features in 
the direction normal to the junction are hidden due to a 
cancellation with the electron velocity, and that the DOS 
in the plane of the junction does not show up due to focal- 
ization effects These two mechanisms cannot explain 
the absence of VHS in c-axis STM/HTS tunnel junctions: 
the HTS materials being quasi two-dimensional have vir- 
tually no dispersion in the tunneling direction, and the 
STM junction, owing to its microscopic size, is qualita- 
tively different from a planar junction and is character- 
ized by a specific matrix element which may not lead to 
focalization effectSi ^^'^^ Besides, the ability of the STM 
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FIG. 1: Typical STM conductance of Bi2223 (Tc = 111 K) 
at T = 2 K (dots). The data is an average of several spectra 
taken at different positions on the same sample, and having 
the same peak-to-peak gap Ap = 38 meV. Error bars give the 
standard deviation of this average. The inset shows the spec- 
trum on a larger energy scale. Also shown are three model 
predictions (see text): free-electrons d-wave BCS (shaded 
area), d-wave BCS with realistic dispersion including VHS 
(dashed line), and d-wave BCS including VHS and coupling 
to the collective mode (full line) . The total spectral weight in 
the energy range of the figure is the same for all curves. 



to probe DOS singularities was recently demonstrated in 
carbon nanotubesi^ 

The solution to this puzzle lies in the coupling to col- 
lective modes. Apart from inducing the dip featurcr^ii^ 
this coupling was also shown to effectively suppress the 
VHS peak in STM spectra. Here we demonstrate that 
the interplay of the VHS and the collective mode con- 
cludes the picture, providing a complete explanation of 
both the missing VHS and the pronounced electron-hole 
asymmetries. Our conclusions are based on a brute force 
fit of the STM spectra and a careful analysis of the model. 
The advantage of this method is to allow for unambigu- 
ous determinations of relevant physical parameters from 
the raw tunneling data, and to provide intuition about 
the relationship between trends in the spectral features 
and parameter variations. 

As a template material we studied the three-layer com- 
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pound Bi2Sr2Ca2Cu308+5 (Bi2223), which has the high- 
est Tc in the Bi-based family. Reproducible spectroscopy 
was recently reported in single crystals of Bi2223 cleaved 
in situ at room temperature, ^-'^ The tunneling conduc- 
tance measured in UHV on an optimally-doped sample is 
shown in Fig.[TJ The spectrum presents the characteristic 
V shape of d-wave superconductors, strong and asym- 
metric coherence peaks, and an asymmetric dip-hump 
structure. We have deliberately selected an optimally- 
doped sample displaying a flat background conductance 
up to high energies (inset in Fig. [T|). For undcrdoped 
samples the spectra acquire an asymmetric background, 
attributed to strong correlation effect Sf^^i^^i^i which can 
conceal other intrinsic asymmetries. By moving to opti- 
mal doping where this background is absent, we can thus 
exclude that the observed asymmetries result from this 
type of correlations. 

For illustrative purposes, we plot in Fig. [T] the predic- 
tion of a conventional free-electrons BCS d-wave model 
(curve defined by the shaded area) i^^i^^ This model fits 
the experimental data well at low energy (< Ap/2), but 
fails to account for the various features present at higher 
energy, in particular the asymmetry of the coherence 
peak height. A much better description of the coherence 
peak height, width, and asymmetry can be achieved by 
taking into account the actual band structure. At en- 
ergies below ~ 200 meV a one-band model turns out 
to be sufficient. We consider the two-dimensional lat- 
tice model ^fc = 2^1 (cos + cos ky) -\- At2 cos kx cos ky + 
2t^{cos2kx + cos2A:y) — /i where U is the i-th neighbor 
hopping energy. For this dispersion the VHS lies at en- 
ergy Cm = —4(^2 — is) — M corresponding to the saddle 
point fc = (tt, 0) = M. We determined the parameters 
of the band through an unconstrained least-squares fit 
of the whole spectrum in the inset of Fig. [1] leading to 
ti = -882, t2 = 239, h = -14, and = -26 meV, 
as well as a d-wave gap Aq = 34.1 meV. It is very en- 
couraging that these numbers determined from nothing 
else than the STM tunneling conductance lead to a Fermi 
surface in semi-quantitative agreement with the one mea- 
sured by photoemission4 In particular, the VHS lies be- 
low the Fermi energy as expected for a hole-doped mate- 
rial. The resulting theoretical curve (Fig. [1] dashed line) 
is very similar to the free-electron model at subgap ener- 
gies, but performs much better up to an energy slightly 
above the coherence peaks. The main effect of the VHS 
is to provide additional spectral weight below the Fermi 
level and thus increase the height of the coherence peaks 
at negative bias. Note that in the absence of broadening 
factors^ the VHS would be visible as a secondary peak 
flanking the main coherence peak (see Fig. [2] below) . 

The "BCS plus VHS" model is nevertheless not sat- 
isfactory above eV ^ 2Ap, where it fails to reproduce 
the significant transfer of spectral weight from the dip 
to the hump, which is strongest at negative bias in the 
experimental spectrum. Generically, such transfers sig- 
nal a strong coupling of the quasiparticles with a col- 
lective excitation, which leads to enhanced damping of 



the former in a limited energy range, and to a simul- 
taneous renormalization of the dispersion. In conven- 
tional superconductors, the electron-phonon coupling is 
known to induce similar features, albeit much less pro- 
nounced, at biases related to the phonon frequencies 
A phonon-based interpretation of the dip-hump in HTS 
has been revivified recently Another candi- 
date is the famous (tt, tt) magnetic excitation^ '^^i^^ com- 
monly known as the "41 meV resonance" — although its 
energy changes from one material to the other and also 
with varying doping. ''■^ Coupling the quasiparticles to this 
collective mode yields a change of the electron self-energy 
which can be expressed to leading order in terms of the 
spin susceptibility Xs(g, w)ii^ Using a parametrization of 
Xs as measured by inelastic neutron scattering, Hoogen- 
boom et al. showed that this model provides a very 
good description of the STM spectra of Bi2212 at several 
dopingSi^*^ Apart from the band-structure parameters ti, 
fi and the d-wave gap Aq , this model has 3 more param- 
eters, namely the resonance energy fig, a characteristic 
length Cs ~ 2a which describes the spread of the collec- 
tive mode around q — {TT,n)r^ and a coupling constant 
9- 

In order to estimate these parameters we again per- 
formed a least-squares fit of the whole spectrum in the 
inset of Fig. [U however keeping the t^'s fixed to their 
values determined previously. This procedure yields 
Aq = 33.9 meV and = —42.4 meV, as well as 
fls = 34.4 meV, in reasonable agreement with the prop- 
erties of the magnetic resonance measured in Bi2223^^ 
The resulting theoretical spectrum matches our exper- 
imental data with amazing accuracy (Fig. [U full line). 
In particular, the model reproduces all of the asymme- 
tries found experimentally between positive and negative 
biases. We would like to stress that these asymmetries 
cannot be understood in models which neglect the band 
structure, as, e.g., in Refs [3ll andfssl. The shape of the 
dip minimum in the theory differs somewhat from exper- 
iment: we shall come back to this below. Fits of similar 
quality have been obtained for many different spectra 
with gaps varying from Ap = 36 to 54 meV. 

The precise interpretation of the theoretical curve in 
Fig. [1] seems complicated due to the interplay of three 
similar energy scales: the d-wave gap Aq, the VHS en- 
ergy Cm, and the collective mode energy Sis, all in the 
30-40 meV range. Still, based on a careful study of the 
model we can identify the origin of each structure in the 
spectrum, as illustrated in Fig. [21 The bare BCS DOS 
Nq(lu) exhibits 5 singularities, namely (a) the V at zero 
energy resulting from the d-wave gap; (6) and (6') the co- 
herence peaks at negative and positive energies {—oJb and 
LUhi respectively); (c) the VHS at energy —lOc below the 
coherence peak, and (c') the weak echo of the VHS at en- 
ergy uic' due to the BCS electron-hole mixing. The inter- 
action with a collective mode leads to inelastic processes 
in which a quasiparticle of momentum k and energy lu is 
scattered to a state with momentum k — q and energy 
LU — Q through emission of a collective excitation with 
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FIG. 2: High-resolution calculation of the DOS in the absence 
[No{u!), dashed line, shifted vertically] and in the presence 
[N{u)), full line] of the coupling to the (tt, tt) resonance. The 
energy and shape of the structures induced by this coupling 
(arrows) is an "inverted image" of the singularities present 
in the BCS DOS (a, b, 6', c, c'), shifted by the mode energy 
Qs- The parameters are as in Fig. (TJ except for Aq — 40, 
= —40, and ila = 30 meV. The inset shows the self-energy 
diagram with the full line representing the BCS Green's func- 
tion and the wavy line the spin susceptibility. 
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FIG. 3: Evolution of the theoretical tunneling conductance 
upon varying Qs (a) or (b). The other model parameters 
are as in Fig. [T] Energies are measured relative to the peak 
maximum at — Ap. The values of Qs and span a larger 
range than their possible variation in the experiment (Fig.[4|, 
but were chosen so with the aim of amplifying the trends. 



quantum numbers (q, 17). The corresponding self-energy 
diagram is sketched in Fig. [21 If the only excitation avail- 
able is a sharp- in-energy mode, all singularities of Nq{uj) 
are mirrored in the self-energy, and exactly shifted by the 
mode energy J7s. Hence the DOS N{uj) including the in- 
teraction with the mode displays 3 pairs of singularities 
indicated by arrows in Fig. [21 the onsets aX uj — ifis, 
below which the quasiparticles do not have enough en- 
ergy to excite a collective mode, a first minimum in the 
dip at —iOb — (resp. LUb' + ^s) corresponding to the 
negative-energy (positive-energy) coherence peak, and a 
second minimum in the dip — echoing the VHS peak in 
No{uj) — which is more pronounced for occupied states 
at —LVc — S^s, but also visible at coc' + ^s- Therefore 
the asymmetry of the dip structure between positive and 
negative biases receives a natural explanation in terms 
of the asymmetry of the underlying BCS DOS, which in 
turn is due to the VHS. The appearance of a double min- 
imum in the dip is a direct consequence of the BCS DOS 
having both a coherence peak at —uJb and a VHS peak 
at —LUc- Such a double minimum is not observed in the 
experimental spectrum of Fig. [Tl At positive bias, the 
various broadening effects^^ are sufficient to smear out 
the two minima into one. On the other hand, we have 
found that if the collective mode has a finite inverse life- 
time of only ~ 6 meV,— then the two minima in the dip 
fade away resulting in a smooth dip also at negative bias 
as observed experimentally. 

The exact relationship between the position of the var- 
ious structures in N{uj) and the parameters of the model 
is not straightforward. The peak maximum, which we 
have denoted by Ap, is not related in any simple way to 
the structures in Nq{uj), although it is numerically close 



to cob- Furthermore uJb is not given by Aq, but by the 
value of the gap function A(fc) = ^(cosfc^; — cosfcj,) 
at the Fermi crossing near the M point. Hence ujb is 
slightly smaller than Aq. The case of lUc is simpler, and 
it can be shown that ojc = y^M+^o- 1^ follows that 
the first minimum in the dip at negative energy lies to 
a good approximation at — Ap — f2s, and the second at 

— ^/^M + ~ ^s- In Fig. [3l we illustrate these two de- 
pendencies by varying fig and independently in the 
model. Our starting point is the spectrum of Fig. [1] re- 
produced in bold in Fig. [3l In order to facilitate the 
comparison we have positioned the spectra relative to the 
negative-energy coherence peak. Varying Qg while keep- 
ing fixed we clearly see that the main change in the 
spectrum is a displacement of the dip and hump relative 
to the peak, consistently with the interpretation given in 
Fig. [21 In particular, the width of the dip at negative bias 
does not depend on ilg. As fJs increases, we also observe 
that the coherence peaks become taller and thinner, while 
developing a shoulder. This shoulder carries part of the 
spectral weight expelled from the dip, and progressively 
exits the coherence peak as the difference in the energy 
scales ils and Ag increases. Fig. [2Ia) further shows that 
ris has not much influence on the electron-hole asymme- 
try of the spectra. In contrast, changing the position of 
the VHS by varying dramatically affects this asymme- 
try. At the lowest considered the spectrum is almost 
symmetric. As the VHS moves toward negative energy, 
the dip at < gets wider (the first minimum in the dip 
does not move, as expected) and the dip at w > dies 
out. Inspection of Fig. [31 also shows that the maximum 
of the hump feature tracks the second minimum in the 
dip, and thus depends on both and ^m- Furthermore 
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FIG. 4; STM conductance spectra of Bi2223 (T^ = 111 K) at 
T = 2 K. Each curve is an average of several spectra taken 
at different locations on the same sample, all having the in- 
dicated peak-to-peak gap Ap. fidip is the energy difference 
between the dip minimum (dot) and the peak maximum at 
negative bias, relative to which voltages are measured. 

the hump gets flattened as |^m1 increases. 

Having clarified the different roles of Vis and in the 
shape of the model DOS, we now come back to experi- 
ment. In Fig. m we plot a series of tunneling conductance 
spectra with peak-to-peak gaps ranging from Ap — 36 
to 54 meV. As a function of increasing Ap, we first ob- 
serve that (i) the dip for occupied states gets wider, (ii) 
the dip for empty states gets weaker, and (iii) the hump 
at negative energy flattens out. These three trends are 
all observed in Fig. Efb), and therefore strongly suggest 
that the increase of Ap goes together with a displace- 



ment of the VHS towards negative energies. This is also 
fully consistent with the idea that local increases of Ap 
in inhomogeneous samples reflect local decreases in the 
hole concentration. Another obvious trend of the data 
in Fig. |4] is that the coherence peaks are reduced with 
increasing gap. As seen in Fig. [^b), this is also consis- 
tent with a shift of the VHS to lower energy. However, 
a look at Fig. El^a) shows that this trend can also be 
ascribed to a decrease in the value of fig- Our calcula- 
tions indeed confirm that Og, as determined by fits to 
the spectra in Fig. [H decreases from 34 to 24 meV with 
increasing Ap. The energy difference between the co- 
herence peak and the dip minimum in the experimental 
spectra, f^dip, also decreases with increasing Ap, but less 
than Vis (f^rom 39 to 35 meV), as can be seen in Fig. H) 
From these numbers it appears clearly that Odip overesti- 
mates Vis by 5 to 10 meV. Recently the energy difference 
between Ap and the inflection point between the dip and 
the hump (extremum in the (Pi /dV^ spectrum) was used 
as an estimate of Vlg in Bi2212,'^- resulting in an average 
value of 52 meV. This same estimate would give a Ap- 
independent result of ~ 57 meV for the data in Fig. 21 
almost a factor of two larger that Vis- 

In summary, we have shown that the van Hove singu- 
larity plays a crucial role in shaping the spectral features 
induced in the STM spectra by the interaction of quasi- 
particlcs with a collective mode. As a result, determining 
the frequency of the mode from STM data is complicated, 
and cannot be done directly from simple structures in the 
dl/dV or d^I/dV^ spectra. This conclusion, obtained 
for a coupling to the (tt, tt) spin resonance, also applies 
to phonon models. 
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